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Abstract. A truncated moment sequence (tms) in n variables and of degree 
d is a finite sequence y = (y a ) indexed by nonnegative integer vectors a := 
(«i, . . . , ct n ) such that «! + ■■■ + ct n < d. Let K C R n be a semialgebraic 
set. The truncated i^-moment problem (TKMP) is: how to check if a tms y 
admits a ii"-measure fi (a nonnegative Borel measure supported in K) such 
that y a = ff^x" 1 ■ ■ -x^dp, for every a? This paper proposes a semidcfinitc 
programming (SDP) approach for solving TKMP. When K is compact, we get 
the following results: whether a tms admits a i<"-measure or not can be checked 
via solving a sequence of SDP problems; when y admits no K"-measure, a 
certificate for the nonexistence can be found; when y admits one, a representing 
measure for y can be obtained from solving the SDP problems under some 
necessary and some sufficient conditions. Moreover, we also propose a practical 
SDP method for finding flat extensions, which in our numerical experiments 
always found a finitely atomic representing measure when it exists. 



1. Introduction 

A truncated moment sequence (tms) in n variables and of degree d is a finite 
sequence y — (y a ) indexed by nonnegative integer vectors a :— (ot\, . . . , a„) 6 N" 
with |a| :=«! + ••• + a n < d. Let K be a semialgebraic set defined as 

(1.1) K={xeR n : gi (x) >0, •■• ,g m (x) >0 }, 

with gi, . . . , g m being polynomials in x. We say a tms y admits a K- measure (a 
nonnegative Borel measure supported in K) if there exists a i^-measure such that 



(1.2) y a = / x a d[i, VaeN" : \a\ < d. 

Jk 

(Here, denote x a := x" 1 ■ ■ ■ x™ n for x = (x\, . . . , x n ) and a = (ax, . . . , a n ).) If 
(|1.2[) holds, we say \x is a representing measure for y and y admits the measure /x. 
The truncated K -moment problem (TKMP) is: How to check if a tms y admit a 
ii'-measure? If it admits one, how to get a representing measure? When K = R n , 
TKMP is referred to as the truncated moment problem (TMP). Let 

(1.3) ^ d :^{y=(y a ):aeN n ,\a\<d}, 

and denote by meas(y, K) the set of all if-measures admitted by y. Let 

(1.4) M d {K) := {y G M n4 : meas(y,K) ^ 0}. 

A measure with finite support is called a finitely atomic measure. A measure /i is 
called r-atomic if |supp(/i)| = r. 
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1.1. Background. Bayer and Teichmann [T] proved an important result: a tms 
y G ^ n ,d admits a K -measure \x if and only if it admits a K -measure v with 
\supp(v)\ < LetR[x] := R[xi, . . . , x n ] be the ring of polynomials in (xi, x n ) 

with real coefficients, and R[x]<j be the set of polynomials in R[x] whose degrees arc 
at most d. Every tms y G ^ n ,d defines a Riesz functional S£ v acting on R[x]<z as 

(1.5) if v ( p Q x Q ) := ^2 PaVa- 

\a\<d \a\<d 

For convenience, sometimes we also denote (p,y) := J£ y (jp). Let Pd{K) denote the 
set of all polynomials in R[x]d that are nonnegative on A. A necessary condition 
for y to belong to & d {K) is that Jzf a is A -positive, i.e., 

JS?„(p)>0 VpePd(K). 

This is because 

^v(p)= / pdti>0 Vpe P d (K),V /j,€meas(y,K). 
Jk 

Indeed, when A is compact, ££ v being A-positive is also sufficient for y to belong to 
Md{K). This was implied by the proof of Tchakaloff's Theorem [30]. A condition 
stronger than Jzfj, being A"-positive is J*f y being strictly K -positive, i.e., 

J2f„(p) > Vp G P d (K), p ^ on A. 

As will be shown in Lemma l2.2[ when A has nonempty interior, a tms y belongs to 
the interior of Md{K) if and only if its Riesz functional Jzf a is strictly A-positive. 
Typically, it is quite difficult to check whether Jzf y is A-positive or strictly K- 
positive. 

A weaker but more easily checkable condition is that the localizing matrix of a 
tms is positive semidefinite if it admits a A-measure. If a symmetric matrix X is 
positive semidefinite (resp., definite), we denote X > Q (resp., X y 0). For a tms 
y of degree 2k, let M k {y) be the symmetric matrix linear in y such that 

^y{p 2 ) = P T M k {y)p Vp € R[x] fe . 

(For convenience of notation, we still use p to denote the vector of coefficients of 
p(x) in the graded lexicographical ordering.) The matrix Mk(y) is called a k-th 
order moment matrix. For h G R[x]2fc, define the new tms h * y such that 

S£ v {hq) = Jfh* v (q) G R[x] 2 fc-do g (? i )- 
The tms h* y is known as a shifting of y. If deg(hp 2 ) < 2k, then it holds that 

(1.6) ^ y {hp 2 ) = p T (M fe _|- dcg( , i)/2 -| (h * y)) p. 

The matrix M%_ rdeg(fe)/21 (h * J/) is called a fc-£/i order localizing matrix of ft. 
If a tms j/ belongs to 3$2k{K), then for i = 0, 1, . . . , m (denote go := 1) 

(1.7) M k ~d t (gi *y)h0 where di = [deg(g>i)/2] . 
This is because for every polynomial p with deg(gip 2 ) < 2k we have 

p T M k _ di (gi * y)p = £fy(giP 2 ) = / g l p 2 dii > 

Jk 
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for all jj, € meas(y, K). Thus, (|1.7[) is a necessary condition for y to belong to 
&2k(K). In general, f|l .7|) is not sufficient for y to belong to 3$2k(K)- However, in 
addition to (|1.7p . if y is also flat, i.e., y satisfies the rank condition 

(1.8) rank M k - dg (y) = rank M k (y), 

then y belongs to Z$2k{K). The above integer d g is defined as 

(1.9) d g := max{l,rdeg( 5i )/2l}. 

l<i<m 

The next important result is due to Curto and Fialkow. 

Theorem 1.1 (|S]). Let K be defined in 11. 1^ and d be even. If a tms y £ ^ n ,d 
satisfies \1. 7| ) and the flat extension condition il.8\) for k — d/2, then y admits a 
unique rank M k {y)- atomic K -measure. 

A nice exposition and proof of Theorem 11.11 can be found in Laurent [30] . Flat 
extensions are not only used in solving TKMPs, but also frequently used for solving 
polynomial optimization problems (cf. [HI HOI 1201 I2"T] ). 

There are other necessary or sufficient conditions for a tms y to belong to &d(K), 
like recursively generated relations. Most of them are about the case of M k (y) being 
singular. We refer to gj [9j [TTJ Qj] for the work in this area. When M k (y) 
is positive definite, there is little work on TKMPs, except for the cases n = 1 (cf. 
0]), or d = 2 or (n,d) = (2,4) (cf. [13]). 

In solving TKMPs, flat extensions play an important role. It is interesting to 
know when and how a tms y € ^C n ,d is extendable to a flat tms w € ^ n ,2k with 
2k > d. We say that y extends to w (or w is an extension of y) if y a = w a for all 
| cat] < d. There is an important result due to Curto and Fialkow [8 about this. 

Theorem 1.2 ([8 ). Let K be defined in Then a tms y € ^ n ,d admits a 

K -measure if and only if it is extendable to a flat tms w € -^ n ,2k with 2k > d and 
Mk-di {g%*w) h for i = 0, 1, . . . , m. 

In view of the above result, the following questions arise naturally: i) If a tms 
y admits no if-measure, how do we get a certificate for the nonexistence? ii) If y 
admits a X-measure, how do we get a representing measure? iii) Preferably, if a 
tms admits a if-measure, how can we get a finitely atomic representing measure? 

1.2. Contributions. This paper focuses on the questions above. We propose a 
general semidefinite programming (SDP) approach for solving TKMPs. When K 
is compact as in (|1.1[) . we have the following results: 

(1) Whether a tms admits a if- measure or not can be checked via solving a 
sequence of semidefinite programs {(SDP)k} (see (|3.2|) or (|3.5[1 ). 

(a) A tms y admits a ^-measure if and only if the optimal value of (SDP)k 
is nonnegative for all k. Consequently, when y admits no if-measure, 
the optimal value of (SDP)k will be negative for some k, which gives 
a certificate for the nonexistence of a representing measure. See The- 
orems 13.21 and 13.31 

(b) When y admits a X-measure, we show how to construct such a measure 
for y via solving (SDP)k for a certain k and using the flat extension 
condition. This works under some necessary and sufficient conditions 
(they are not far away from each other) . See Theorems 13.21 and 13.61 
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(2) We propose a practical SDP method for finding flat extensions, and thus 
provide a way for constructing finitely atomic representing measures. It is 
based on optimizing linear functionals (see (|4.1[l ) on moment matrices and 
consists of a sequence of SDP problems. In our computational experiences, 
the method always produced a flat extension if it exists. We derive a bit 
of supporting theory for this fact. When a tms admits a _ftT-measure, we 
prove that for a dense subset of linear functionals, the method asymptoti- 
cally produces a flat extension. When a tms does not admit a -ftT-measure, 
we prove that this sequence of SDP problems will become infeasible after 
some steps. This method is also applicable when K is noncompact. See 
Subsection 4.1 and Theorem 14.31 

The results described in (1) are given in Section[3l and those in (2) are in Section[4] 
Section 2 presents some background for proving these results. 

2. Preliminaries 

2.1. Notation. The symbol N (rcsp., K.) denotes the set of nonnegative integers 
(resp., real numbers), and R" denotes the nonnegative orthant of R n . For t € K, 
\t] (resp., [t\) denotes the smallest integer not smaller (resp., the largest integer 
not greater) than t. The [x]d denotes the column vector of all monomials with 
degrees not greater than d: 

[x] d = [I x\ ■■■ x n x\ xix 2 xf x^~ 1 x 2 xtV ■ 

For a set S C W 1 , \S\ denotes its cardinality. The symbol int(-) denotes the interior 
of a set. For a matrix A, A T denotes its transpose; if A is symmetric, X m i n (X) 
denotes its minimum eigenvalue. For u e R N , \\u\\ 2 := V u T u denotes the standard 
Euclidean norm, and B(u,r) := {x £ l n : ||x — u|| 2 < r} denotes the closed ball 
with center u and radius r. The • denotes the standard Frobenius inner product 
in matrix spaces. For a matrix A, \\A\\p denotes the Frobinus norm of A, i.e., 
\\A\\p — \/Trace(A T A) . The zero set of a polynomial q is denoted by Z(q). For 
a tms w, w\t denotes the subsequence of uu whose indices have degrees not greater 
than t, i.e., w\t is a truncation of w with degree t. For a measure /1, supp(/i) denotes 
its support. A polynomial / 6 M.[x] is said to be sum of squares (SOS) if there exist 

fi,...,fk & KM such that / = ff H h /|. The set of all SOS polynomials in n 

variables and of degree d is denoted by S„ ^. The symbol ^{K) denotes the space 
of all functions that are continuous on a compact set K, and || • denotes its 
standard cx>norm. 

2.2. Truncated moments. For a compact set K, a tms admits a measure sup- 
ported in K if and only if its Riesz functional is if-positive. This result is implied 
in the proof of Tchakaloff's Theorem 30 (cf. [13 ). 

Theorem 2.1 (Tchakaloff). Let K be a compact set in W l . A tms y admits a 
K -measure if and only if its Riesz functional Jzf y is K -positive. 

In the following, we characterize the interior of &d(K) via strict if-positivity of 
Riesz functionals. 

Lemma 2.2. Let K C W 1 be a set with int(K) ^ and y be a tms in ^ n ,d- Then 
y belongs to int(3%d{K)) if and only if S£ y is strictly K -positive. 
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Proof. " => " Suppose y belongs to int {3%d{K)). Let £ £ ^ n ,d be the tms generated 
by the standard Gaussian measure restricted to K. Then, for every p £ Pd{K) with 
p\k ^ 0, we must have > because int(.K') 7^ 0. If e > is small enough, 

the tms v := y — e£ belongs to Md{K) and it holds that 

jSf„(p) = JSf„(p) + cJSf c (p) > e^ c (p) > 0. 

This means that Jz? y is strictly X-positive. 

" •£= " Suppose ££ y is strictly X-positive. Lemma 2.3 of [13] implies that, for 
some 5 > 0, the is strictly if-positive for all w £ M n ,d with — y\\ 2 < S. 
Then, by Theorem 2.4 of [T3], every such a tms w, including y, belongs to &d(K). 
This implies that y belongs to int (Md{K)). □ 

A tms w £ ^ n ,2k generates the moment matrix M k {w). Recall that, for a 
polynomial /, we still denote by / the vector of its coefficients. The vector / is 
indexed by exponents of x a . The matrix vector product M k (w)f is defined in the 
usual way, i.e., 

M k (w)f = & w (f[x] k ). 
We say p £ kerM k (w) if M k (w)p = 0, i.e., Jz? w (px a ) = for every |a| < k. 

An ideal of R[x] is a subset / C R[x] such that I + I C I and p ■ q £ I for 
every p £ I and q £ R[x]. Given p\, . . . ,p m £ R[x], denote by (pi, ■ • • ,p m ) the 
ideal generated by p%,. . . ,p rn . If I = (pi, . . . ,f> m ) and every ^ £ ker M k (w), then 
/ C (kerMfc(w)), the ideal generated by polynomials in kerMfe(w). 

Lemma 2.3. Let w £ ./#„.2fc, h,p £ R[x] 6e smc/i i/iai -ff := Affe-rd e g(h)/2l {h*w) > 
0. Ifp,q are polynomials with p £ ker H and 

deg(pq) <k- \deg(h)/2] - 1, 

i/ien we have pq £ ker H. 

Proof. Let z = h * w. Then z is a tms in ^ n ,2k-dcg(h) an d the moment matrix 
A jr fe-[dog(/i)/2] ( z ) b 0. The conclusion is implied by Lemma 3.5 of [TB] (also see 
Lemma 5.7 of [?T] or Theorem 7.5 of [5]). □ 

2.3. Quadratic module, preordering and semidefinite programming. For 

the semialgebraic set K defined in (jl.ll) . its fc-th truncated quadratic module Q k (K) 
and truncated preordering Pr k (K) are respectively defined as 



(2.1) Qk(K):=lJ2^ 

{ i=0 



each deg(<7i(?i) < 2k 
and o~i is SOS 



deg{a u g u ) < 2k 
each a v is SOS 



(2.2) Pr k {K) := I £ a v9u 

{ ^e{o,i}™ 

i.9v '■= 9 1 1 ' ' ' 9 V m ■) The quadratic module and preordering of K are then defined 
respectively as 

Q(K) = |J Q k (K), Pr(K) = (J Pr k (K). 

k>0 k>0 

The definitions of Q k (K) and Pr k (K) depend on the set of defining polynomials 
<7i, . . . ,g m , which are not unique for K. Throughout the paper, when Q k (K) or 
Pr k {K) is used, we assume that gx,...,g m are clear in the context. In (jl.l|) . if K 
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is defined by using polynomial equalities, like h(x) — 0, then it can be replaced by 
two inequalities h(x) > and —h(x) > 0. 

Theorem 2.4. Let K be as in and p E M.[x] be strictly positive on K. 

(i) (Schmiidgen, |26] ) If K is compact, then we have p E Pr(K). 

(ii) (Putinar, |24] ) If the archimedean condition holds for K (a set {x : q(x) > 
0} is compact for some q E Q{K) ), then we have p E Q{K). 

Theorem 2.5 (Real Nullstellensatz). Let K be defined in iTP]) . If f E R[x] 
vanishes identically on K , then —f 2i E Pr(K) for some integer I > 1. 



Remark: Theorem 12.51 is a special case of the so-called Positivstellensatz [27] , and 
the set K there does not need to be compact. 
Let I(K) be the vanishing ideal of K, i.e., 

I{K) := {h E R[x] : h(u) = V u E K}. 

By Theorem^ if / E I(K), then -f 2i E Pr k (K) for some k,L This fact will be 
used in our proofs later. 

The sets Qk(K) and Prk(K) are convex cones. Their dual cones lie in the space 
~^n,2k- The dual cone of Prk(K) is defined as 

Pr k (K)* = {yE Jt n . 2k : (p, y) > Vp e Pr k (K)}. 

The dual cone Q k (K)* is defined similarly. It is known that (cf. [161 119) ) 

Qk(K)* ={i/6 ^n,2k ■ Mh-di(9i *y)h0 for i = 0, 1, . . . , m} , 
Pr k {K)* = {yE ^ 2k : M k ^ dv (g v * y) h E {0, 1}™} . 

(c?, y := |~deg (g^)/ 2] .) Given tms' c, ai, . . . , at E ^ n ,d and scalars 6i, . . . , b t , consider 
the linear conic optimization problem 



(2.3) 



(P, c) 

(p,Oi) = b» (1 < i < t), P ER[x} d nPr k (K). 
Its dual optimization problem is 




(2.4) 



max &1A1 + • • • + btXt 

s.t. w\d = c — Aiai \ta t , 

M k -d„(g»*w)h0Vv G {0,l} r 



The optimization problems (|2.3p and (|2.4p are reducible to SDP problems (cf. [TBI 
[HI [21]). Any objective value of a feasible solution of (12 .3[) (resp., (|2.4[l ) is an upper 
bound (resp., lower bound) for the optimal value of the other one (this is called 
weak duality). If one of them has an interior point (for (|2.3I) it means that there is a 
feasible p lying in the interior of R[x} d nPr k (K), and for (12.41) it means that there is a 
feasible w satisfying every M k -d„ {gv * w) >- 0), then the other one has an optimizer 
and they have the same optimal value (this is called strong duality). Similar is 
true if Pr k (K) is replaced by Q k (K). We refer to [TBI Ull [21] for properties of 
SDPs arising from moment problems and polynomial optimization. In [17] Lasserre 
proposed a semidefinite programming approach for solving the generalized problem 
of moments. We refer to [31] for more about SDP. 
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3. Checking Existence of Representing Measures 

This section discusses TKMPs when if is a compact semialgebraic set defined in 
(|1.1[) . It gives a semidefinite approach for checking if a tms admits a representing 
measure and proves its properties alluded to in £jl.2f 1). 

3.1. A certificate via semidefinite programming. Our semidefinite approach 
for solving TKMPs exploits the following basic fact. 

Proposition 3.1. Let y be a tms in ^ n ^ d and K be a compact set defined by U.l\) . 
Then y admits no K -measure if and only if there exists a polynomial p such that 

(p,y)<0, P eR[x] d nPr(K). 

Under the archimedean condition for K, the above is also true if Pr(K) is replaced 
byQ(K). 

Proof, (i) The "if" direction is obvious. It suffices to prove the "only if" direction. 
Suppose y does not belong to & d (K)- Then, by Theorem 12. 11 the Riesz functional 
J£ y must achieve a negative value on Pd(K ), say p £ Pd(K), such that ^f y (p) < 0. 
Then, for a small enough e > 0, the polynomial p := p + e also satisfies J£ y {p) < 0. 
Since p is strictly positive on the compact set K, by Theorem I2.4( i). we have 
p 6 Pr(K). This p satisfies (p,y) < 0. The proof is same when the archimedean 
condition holds (applying Theorem l2.4f ii)). □ 

In the following, we show how to apply Proposition 13. II to check whether a tms 
y belongs to £% d (K) or not. Choose a tms £ € M d {K) such that _Sfj is strictly 
if-positive. For an integer k > d/2, consider the optimization problem 

(3.1) \ k := min (p,y) s.t. (p,£) = 1, p £ R[x] d H Pr k {K). 

v 

Its dual optimization problem is {d v :— |~deg(<7„)/2]) 

!Afe := max A 
s.t. w\d=V-*£, 
M k ^ dij (g„*w)hOVv£{0,l} m . 

Both the primal (|3.1[) and dual Q3.2p are reducible to SDP problems, and they 
are parameterized by an order k. They can be solved efficiently, e.g., by software 
GloptiPoly [15| and SeDuMi [28] . 

A nonnegative Borel measure [i is called (K,d)- semialgebraic if supp(/x) C K 
and there exists q £ M.[x] d fl Pr(K) such that 

supp(^i) C {u £ K : q{u) = 0}, q ^ on K . 

If a measure [x is not (K, (i)-semialgebraic, then 

suppOu) C K n Z(q), qeR[x] d n Pr{K) q = 0onK. 

Note that not every if-measure is (K, d)-semialgebraic. For instance, for K being 
the unit ball, the probability measure uniformly distributed on K is not (K, d)- 
semialgebraic, because every polynomial vanishing on it must be identically zero. 
So, being (K, d)-semialgebraic is a bit restrictive for a measure ^i, as it implies 
supp(/i) has Lebesgue measure zero when K has nonempty interior. 
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Theorem 3.2. Assume the set K defined in fZHP is compact, the tms £ belongs to 
Md{K) and its Riesz functional Jz?£ is strictly K-positive. Let y, Xk, Afe be as above. 
Then the sequence {Afe} is monotonically decreasing. Let Aoo := linife._i.oo Afe € 
KU {— oo}. We also have: 

(i) For k big enough, Xk = Xk and 113.2]) has a maximizer. 

(ii) The tms y belongs to &d(K) if and only if Xk > for all k. If Aoo > — oo. 
then the shifted tms y := y — Aoo£ belongs to 3%d{K)- 

(iii) Assume int(K) ^ 0. Then, Aoo > — oo, and Xk — Aoo for some k if and 
only if there exists a (K, d)-semialgebraic measure fi G meas(y, K). 

(iv) Suppose, for some ko, \3. S\) has an optimizer (Afc ,w*) with Xk > and 
w*\t(t > d) flat. Then y belongs to &d{K) and Xk = Xk for all k > ky. 

Proof. Since the feasible set of (|3.2[) is shrinking as k increases, the sequence of op- 
timal values {Afe} must be monotonically decreasing. So its limit Aoo :— limfe_>.oo Afe 
exists (would possibly be — oo). 

(i) Since K is compact, there exists R > big enough such that, for all a E 
N™ with |a| < d, the polynomial R — 1 + x a is positive on K and belongs to 
Prk(K) for some k, by Theorem I2.4f i). For e > small enough, we also have 
R + c + (1 + a)x a e Pr k (K) whenever |c |, |ci| < e, for all a € N n with \a\ < d. 
This implies that the polynomial p :~ J2\ a \<d(^ + xa ) ues m * ne interior of the 
cone M.[x]d fl Prk(K). Thus, for k big enough, the primal problem (|3.1I) has a 
feasible point lying in the interior of M.[x]d fl Prk(K). Hence, the primal and dual 
optimization problems (|3.1j) and (|3.2[) have the same optimal value, and (|3.2p has 
a maximizer (see the comments following (|2.4p ). 

(ii) If the tms y belongs to £%d{K), then Afe > for all k, because A = is 
feasible in (|3 . 2[) for every k. If y does not belong to &d(K), by Proposition 13.11 
there exists p £ R[x]d n Prk(K) such that (p, y) < 0, for some k. Generally, we can 
assume p ^ on K (otherwise, replace p by p + e for a tiny e > 0). By the strict 
iC-positivity of £, we can normalize p as (p, £) = 1. Hence, Afe < Afe < 0. This 
shows that the first statement of item (ii) is true. 

The second statement of item (ii) is implied by the first one. This is because, in 
(|3.2p . if y is replaced by y, then the resulting optimal values are all nonnegative. 

(iii) By item (i), for big k, Xk = Xk is at least the optimal value of 

(3.3) min (p,y) s.t. (p, £) = l,p G P d (K). 

p 

Since mt(K ) ^ 0, p ^ on K if and only if p is not identically zero. Thus, the 
strict i^-positivity of Jz?£ implies 

(3.4) 1 = > ei||p|| 2 , for some a > 0. 

So, the feasible set of (|3.3[) is compact and has a finite minimum value A > — oo. 
By the monotonicity of {Afe}, we know Afe > A for every k and its limit Aoo is finite. 

"if direction: If there is a (K, ci)-semialgebraic measure fi € meas(y, K), then 
for some k we can find ^ q e R[x]d n Prk{K) such that supp(^i) C K fl -2-(g). 
Normalize g as (<z,£) =1, then 

{q, y) = (q, y) + ^oo = / qd^ + Xoc = x^. 

JK 

So, the polynomial q is feasible for (|3.1|) . and Afe < A^. By weak duality, Afe < 
Afe < Aoo- Since Afe > A^, for all k, we must have Afe = Aoo- 
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"only if direction: If Afc = Aoo for some ko, then Afc = Aoo for all k > ko. Fix 
such a k. Since int(A) ^ 0, the feasible set of (|3.1|) is compact (by (|3.4|l ) and has 
a minimizer p* ^ 0. As shown in (i), if k is big enough, Afc = Afc and (|3.2p has a 
maximizer, say, (\k,w*). Then, for every p, g meas(y, K), it holds that 



= (p*,y-X k & = ip\y) = / P*dfx. 

JK 

Since p* is nonnegative on A, supp(/^) C Z(p*) and /i is (A , d)-semialgebraic. 

(iv) Note that y = w*\d + \k £, and £ admits a A'-measure, say, v. If w*\t(t > d) 
is flat, then w*\t admits a finitely atomic A-measure /Lt, by Theorem 11.11 Since 
Afc > 0, we know [i + Xk v is a representing A-measure for y. So, y belongs to 

• *',/; A i. 

Clearly, for every k, the projection of the feasible set of (|3.2[) into the A-space 
contains the set 

A={A: y-A£g^ d (A)}. 

Thus, every optimal value Afc is greater than or equal to the optimal value 

A* := max A s.t. A e A. 

So, we know Afc g A from the above. Hence, Afc > A* > Afc for all k. By the 
decreasing monotonicity of {Afc}, we know Afc = Afc for all k > kg. □ 

When int(A) ^ 0, as we can see in the proof of item (iii) of Theorem l3.2[ a mea- 
sure in meas(y, A) is (A, d)-semialgebraic if and only if all measures in meas(y, A) 
are (A, (i)-semialgebraic. 

In (|3.2p . all the cross products g v :— g^ 1 ■ ■ ■ g 1 ^ are used, which might be incon- 
venient in applications. So, we consider a simplified version of 



(3.5) 



Afc := max A 

s.t. w\d = y-\£, 

Mk-di (gi * w) h (0 < i < to). 



The dual optimization problem of (|3.5[) is 

(3.6) A^. := min fay) s.t. (p, £) = 1, p g R[x] d n Q fc (A). 

p 

We have the following analogue to Theorem 13.21 

Theorem 3.3. Let A, £, y, y, Afc, Aoo &e same as m Theorem \3J^ and Afc, A^ 6e as 
above. Suppose the archimedean condition holds for A . Then, the sequence {Afc} is 
monotonically decreasing. Let Aoo :— linifc^oo Afc g R U {— oo}. We also have: 

(i) Aor k big enough, Afc = X' k and US. 5\) has a maximizer. 

(ii) The tms y belongs to &d{K) if and only if Xk > /or aZ/ /c. 

(iii) A always holds that Aoo = Aoo- If int(K) ^ 0, then Aoo > — oo, and 
Afc = Aoo for some k if and only if there exists fi g meas(y, A) such that 

supp{n) C A n Z(a) /or some q g n Q(A), q\ K ^ 0. 

(iv) Suppose, for some ko, the optimization problem 113.5]) has an optimizer 
(Afc ,w*) with Afc > and w*\ t (t > d) /Zai. Then, y belongs to Md(K) 
and Afc = Afc /or all k > kg. 
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Remark: By Theorem l3.3l (ii). if y does not belong to Md{K), then Afc < for some 
k, and for A = there is no w G ^ n ,2k satisfying 

(3.7) w\ d = y, M k (w) h 0, M fe _ rfl ( 5l * w) h0, M k _ dm (g m *w) h 0. 
So, we have y g" Md{K) if and only if (|3.7|) is infeasible for some A;. 

Proof of Theorem \3.3[ The decreasing monotonicity of {Afc} holds because the fea- 
sible set of (|3.5[) shrinks as A; increases. 

(i) As in Theorem 13. 21 under the archimedean condition, for k big enough, we 
can similarly show that (|3.6p has a feasible point lying in the interior of the cone 
E.[x]d n Qk{K). Thus, (|3.5[) and (|3.6p have the same optimal value, i.e., Afc = X' k , 
and (|3.5p has a maximizer when fc is sufficiently large. 

(ii) This can be proved in a way similar to item (ii) of Theorem 13.21 under the 
archimedean condition on K. 

(hi) We first show that Aoo = Aoo- Clearly, Afc > A fc for all k. Thus Aoo > A^. 
If Aqo — ~oo, then we are done. If > — oo, we need to show A^, = Aoo- If not, 
seeking a contradiction, we suppose otherwise Aoo > Aoo- Let y := y — Aoo£- If we 
replace y by y in (|3.5j) . then all its optimal values are nonnegative, and thus we 
have y <E &d(K) by item (ii). However, by Theorem 13.21 (ii), we get y ^ £%d(K), a 
contradiction. This is because if we replace y by y in (|3.2[) , then its optimal value 
is Afc — Aoo which is negative for big k. So, we must have Aoo = A^. 

The second statement of item (hi) can be shown in a similar way as for item (iii) 
of Theorem 13.21 The only difference is to apply Proposition 13.11 with Q(K), under 
the archimedean condition on K. 

(iv) The proof is the same as for item (iv) of Theorem 13.21 □ 

Example 3.4. Consider the tms y £ ^#2,6 : 

(28; 0,0; 1.1, 0,3.4; 0,0, 0,0; 1.1, 0,1.2, 0,1.6; 0,0, 0,0, 0,0; 28, 0,3.4, 0,1.6, 0,1.2). 

Its 3rd order moment matrix M^(y) is positive definite. Let K — {x € M 2 : ||x||| < 
25}. Choose £ to be the tms in ^2,6 induced by the probability measure uniformly 
distributed on the unit ball. For k = 3, 4, 5, we solve (|3 . 5|) (which is same as (|3.2[) . 
since K is defined by a single inequality) using GloptiPoly |15j . and get optimal 
values Afc numerically as 

A 3 w 0.3702 > 0, A 4 w 0.0993 > 0, A 5 ~ -0.2370 < 0. 

By Theorem 13.21 or 13.31 we know this tms does not admit a if-measure. □ 

By Theorem l3.2[ when K is compact, a tms y belongs to &d(K) if and only if Afc 
is nonnegative for all fc, which is difficult to check in applications. Similarly, it is 
also difficult to check Afc > for all k. However, the items (iv) of Theorems 13.21 and 
13.31 provide a certificate for checking the membership y 6 £%d{K). This is because 
we can easily check the flat extension condition (|1.8I) as follows. Let (A w*) be an 
optimal pair for (|3.2p or (|3.5p . If a truncation w*\ t (t > d) is flat and Afc > 0, a 
JC-measure representing y can be constructed from the relation y = w*\d + Afc£, 
because the tms £ admits a measure by its choice. If we have [i € meas(w* , K) and 
v e meas(^, K), then /i + X k v belongs to meas(y, K). 



Throughout the paper, the moments of a tms are listed in the graded lexicographical ordering, 
and moments of different degrees are separated by semicolons. 
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Example 3.5. Consider the tms y £ ^2,6- 

(1; 7/50, 1/50; 2/5, 0, 2/5; 91/1250, -6/625, 42/625, 37/1250; 
6973/25000, -42/3125, 1777/25000, 42/3125, 6973/25000; 1267/31250, 
-222/15625, 504/15625, 72/15625, 546/15625, 781/31250; 23713/100000, 
-42/3125, 2929/100000, 0, 2929/100000, 42/3125, 23713/100000). 

Let K be the 2-dimensional unit ball; thus d g = 1. The 3rd order moment matrix 
M 3 (y) is positive definite. Let £ be the tms in ^2,6 induced by the uniform prob- 
ability measure v on the unit ball. We solve (|3.2[) by GloptiPoly [15]. For k = 5 
we get the optimal value A5 « 1.2 • 10~ 7 . The computed optimal w* is flat, because 
rankM 4 (w*) = rankA/ 5 (w*) = 10. By TheoremO we know w* admits a unique 
10-atomic measure /i. □ 



For a fiat tms, a numerical algorithm is given in [14] to find the support of its 
finitely atomic representing measure. In Example 13.51 by using this algorithm, we 
get the support of the 10-atomic measure fj, admitted by w* there: 



From y = w*\e + As£, we know fi + A5 • v is a ball-measure representing the tms y 
there. 

3.2. Finding representing measures via shifting. The sign of Aoo is critical 
in checking the membership in Md{K). When int(-£sT) ^ 0, by Lemma \2. 21 a tms y 
belongs to int (&d(K)) if and only if its Riesz functional Jzf a is strictly if-positive. 
By Theorem 12. 1[ for compact K, we know ^%d(K) is closed. Thus, for compact K 
with 'mt(K) 7^ 0, we can check the membership in £%d{K) as: 

• If Aqo > 0, y lies in the interior of £%d{K)] 

• if Aqo = 0, y lies on the boundary of &d(K); 

• if Aqo < 0, y lies outside .^d{K). 

When Aoo > —00, the shifted tms y :~ y — Aoo£ always admits a if-measure. When 
Aoo > 0, if we have \i <E meas(y, K) and v 6 meas(^, K), then /i + Aqo^ is a K- 
measure representing y, because of the relation y = y + Aoo£- Thus, it is enough to 
investigate how one can get a if-measure representing y. 

When int(K) ^ 0, Theorem I3.2f iii) implies that Aoo > —00 and A& = Aoo 
is achieved at some step k if and only if a measure representing y is (K, d)- 
semialgebraic. We are mostly interested in this case, since it is not possible to 
solve (|3.2[) for infinitely many k's. Therefore, in the following, we assume y admits 
a (K, (i)-semialgebraic measure. 

For an optimizer w* of (|3.2[) . the kernel of the moment matrix Mk(w*) is a 
useful tool in constructing supports of representing measures and their vanishing 
ideals. In particular, the kernel kerM/ c (w*) is of the greatest interest when Mk{w*) 
achieves the maximum rank over the set of all optimizers of Q3.2p . This approach 
was introduced by Lasserre, Laurent and Rostalski [18) when they compute zero- 
dimensional real radical ideals via semidefinite relaxations. For more details, we 
refer to the surveys [21] by Laurent and [22] by Laurent and Rostalski. In the 
following, we use this approach. 

For each fc, denote by 0pt k (K) the set of optimizers of (13.21) . 

2 The ranks here are evaluated numerically. We ignore singular values smaller than 10 — 6 when 
evaluating ranks. The same procedure is applied in computing ranks throughout this paper. 



{±(1,0), ±(0,1), (±1/2, ±1/2), (4/5,-3/5), (3/5,4/5)}. 
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Theorem 3.6. Let K be defined in U.l\) (not necessarily compact), and y, £, y, X k , A c 
be same as in Theorem VS.'A Suppose \ ka = Aco and there exist /i £ meas(y, K) and 
7^ q G M[x]d fl Pr kl (if), /or some ko, k\, such that 

(3.8) supp(n) C F:={neJf: g(w) = 0}. 

Lei fc > max(fc , fci). 1/ Dpt k (K) 7^ 0, t/ien /or eac/i (A, w) G Dpt k (K) we have: 

(i) For k sufficiently large, it holds that I(V) C (ker Mk(w)) . 

(ii) Let r = max{ rankM k (v) : (A, v) G Dpt k (K)}. If rankMk(w) — r, then 
(kevM k (w)) C I(supp(fi)). 

(iii) // I V| < oo, £/ie7i w|2fc-2 *s /?a£ /or sufficiently large. 

(iv) 7/ |swpp(/i)| = 00 and rankM k (w) — r, then w\2t can not be flat for all 
< t < k. 

Proof. Since q is in Pr kl (K), there exist SOS polynomials <j v such that 

Write each a v = J2j ®tj- F° r every k > max(fco, fci), it holds that (note w\ d = y) 




Because M k _ dv {g v * w) h for all u, we have M k _ dv {g v * w)6 v j = for all v,j. 

(i) Let {hi, . . . , h r } be a Grobner basis of the vanishing ideal I(V), under a total 
degree ordering. Then, each hi vanishes on V. By Theorem 12. 51 for each h i} there 
exist an integer I > 1, a polynomial </> and SOS polynomials tp u such that 

(3.9) hf + qcj> + <Pv9» = 0- 

Write <f) — <j)\ — (f>2 for some (f>i,(f>2 G R[x}. From v j G kerM d _ d _(g„ * w), by 
Lemma 12.31 we must have v ^<p s G ker M d - dv {g v * u>) if 

deg(^,^ s ) + rdeg( ff ,)/2l < fc - 1 (s = 1, 2). 

This implies that 

{q<f>,w)= Y (^u,j<t>i) T M k ^ du (g u *w)(d Utj (l)i) 
^e{o,i} m j 

- {6 vd fc) T M k - du (g v * w)(6 Vtj ct> 2 )} = 0. 

Combined with (|3.9I) . this gives 

(M fc H^,^)+ J] =0. 
ye{o,i} m 

Since M k - dv {(j v * w) ^ for all z/, we know ((/?„<?„, u>) > for all i/, which implies 
M k (w)hf = 0. By an induction on we can get hi G ker Mfc(u>) (cf. Lemma 3.9 of 
PU). So, /(F) C (kerM fe (w;)) for k big enough. 

(ii) Let v G .«# n ,2fc be the tms induced by /i, i.e., w = J[x]2ferf/U. Then (Afc,«) 
belongs to Dpt k (K). The rank M k (w) being maximum over 0pt k (K) implies 

rank M k (w) > rank M k ((w + v)/2) = rank M k (w + v). 



A SEMIDEFINITE APPROACH FOR TKMP 



13 



Since both M k (w) h and M k (v) y 0, it holds that 

ker M k (w + v) = ker M k (w) H kerM k (v). 
The relation rank M k (w) > rank M k (w + v) and the above imply 

kerMfc(u>) = ker M k (w + v) C ker M k (v). 
For all / £ kerMfc(w), we have / € ker Mfc(v) and 

= f T M k (v)f= f fdfi. 

Jk 

Thus, / vanishes on supp(/i), and (ker M k (w)} C 7(supp(/^)). 

(hi) When |V| < oo, the quotient space R[x]//(y) is finitely dimensional. Let 
{bi, . . . , b t } be a standard basis for it. Then, every monomial x a can be written as 

x a = r(a) + ^2 fchi, deg((j>ihi) < \a\, r(a) € span{b\, . . . ,b t }. 

If k is big enough, we have hi £ ker M k (w) for all i, by item (i). We also have 
4>ihi £ ker M k (w) for all i if \a\ < k — 1, by Lemma \2. 31 Thus, x a — r(a) belongs 
to kerMfc(w) for all |a| = k — 1. So, if 

k - 1 - d 9 > max{deg(6i), . . . , deg(6 t )} 

is big enough, every a-th (|a| = k — 1) column of M k (w) is a linear combination of 
the /?-columns of M k (w) with \j3\ < k—l — d g . This means that the tms generating 
the moment matrix M k -\(w) is flat, i.e., the truncation w\2 k -2 is flat. 

(iv) For a contradiction, suppose w\2t is flat for some < t < k. Then, for 
every |a| = t, there exists a polynomial ip a e R[x] t _i such that x a — tp a belongs to 
ker M t (w). This implies that 

= (x a - <p a ) T M t (w)(x a - <p a ) = (x a - tp a ) T M k (w)(x a - <p a ). 

Since M k (w) >z 0, we have x a — tp a £ ker M k (w). By item (ii), we have x a — Lp a £ 
I(supp(/i)). Then a simple induction on \a\ shows that 

x a = ip a mod 7(supp(//)) 

for some ip a £ R[x]t_i, whenever \a\ > t. That is, every high degree polynomial is 
equivalent to a polynomial of degree not bigger than t — 1 modulo 7(supp(/i)). So, 
7(supp(/i)) is zero-dimensional, i.e., dimR[x]//(supp(/i)) < oo, which contradicts 
|supp(/i)| = oo (cf. Proposition 2.1 of Sturmfels [H])- □ 

Remarks: a) In Theorem I3.6[ we do not need to assume K is compact, b) The 
condition (|3.8p implies J£y(q) — and the Riesz functional Jzf^ is not strictly K- 
positive. Thus, the tms y lies on the boundary of Md{K). This is an advantage of 
using y instead of y. c) When (|3.8I) holds and |V| < oo, the truncation w*\2 k -2 is 
flat for k big enough, but w* itself might not be flat. 

3.3. Tms with a finite if -variety. In this subsection, we consider tms' whose 
associated algebraic varieties have finite intersection with K. Let do = L^/2J- The 
algebraic variety associated to a tms y £ ^ n .d is defined as (cf. |13| ) 

(3.10) V{y):= f] {x £ R™ : p{x) = 0}. 

pEkci- M do (y) 

Then, we define the K-variety associated to y as 

(3.11) V K (y) := V(y)nK. 
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When the moment matrix Md a {y) is singular (thus the variety V{y) is a proper 
subset ol R™, like a finite set or a curve, etc.), there exists work on discussing 
whether y admits a measure. We refer to [71 EJ H2] ■ 

Clearly, for Vif(y) to be a proper subset of K, the moment matrix Md (y) must 
be singular, and any A > is not feasible in (13.2[) . To see this, suppose otherwise 
a pair (A, w) with A > is feasible in (|3.2[) . When Vk(v) K, there must exist 
p £ ker Md (y) such that p ^ on K . The relation w\d = y — A£ implies 

Md (y) = M do (w)+\M do {0- 

The above moment matrices are all positive semidefinite. Thus, from p £ ker Md (y) 
and A > 0, we know p £ kerMd (£) and ^{p 2 ) = 0. This contradicts the strict 
JT-positivity of £. Therefore, when Vk(v) ^ K, the optimal value Xk of (|3.2I) must 
be zero if y belongs to 3£d{K), and (|3.2j) is equivalent to the problem 

(3.12) find w£Jt ni 2k s.t. w\d = y,M k -d l ,(9v*w)h0Vv£{0,l} m . 

Theorem 3.7. Let K be defined in (not necessarily compact), and y be a tms 

in ^ nt d- Suppose y belongs to 3%d(K) an d \^k(v)\ < oo. If k is sufficiently large, 
then every w satisfying \3.12)) has a flat truncation w\2k— 2- 



Proof. As we have seen in the above, every optimal Xk in (|3.2[) is zero. Thus, 
Aoo = and y = y — Aoo£ = y. Choose a measure /i £ meas{y 1 K) and a basis 
{pi, . . . ,p r } of ker M do (y)- Let q = p\ -\ 1 p 2 r £ R[x] d - From 

/ qd/j, = y2 pl d v = y2pI M d {y)pt = Q, 

Jk *=i Jk *=i 

we know supp(//) C Z(q) n K. The finiteness of Vk{v) implies the set Z{q) n K is 
finite. Then, this theorem follows from item (iii) of Theorem 13.61 □ 



Remark: As one can see in the proof, Theorem l3.7l is an application of Theorcm l3.6l 
If the moment matrix Md {y) is singular and the tms y belongs to &d(K) and 
Vk(v) K, we have already seen that the optimal value of (|3.5|) must be zero, 
and hence (13.5[) is equivalent to the feasibility problem (13.12[) (see the argument 
preceding Theorem 13 .7p . 

Example 3.8. Consider the tms y £ ^#2,4: 

(1; 0, 0; 1, 0, 1; 0, 0, 0, 0; 1, 0, 1, 0, 1). 

The set K = R 2 , do = 2, and the moment matrix Md (y) is singular (it has rank 
4 and size 6x6). The K- variety is defined by x\ = x\ = 1. For k = 4, we solve 
(|3.12[) and get a feasible tms w £ Jtis- Its truncation w\§ is fiat (while w itself is 
not), and it, as well as y, admits a 4-atomic measure supported on {(±1, ±1)}. □ 

3.4. The case K = R". When K = R n is the whole space, we can similarly apply 
the approach described in the preceding subsections. However, since R™ is not 
compact, weaker conclusions can be made. 

Let y be a tms in ^ n ,d and do = \_d/2\ . If y is extendable to a flat tms w £ ^ n ,2k, 
then we can get a finitely atomic representing measure for y, by Theorem ll.il In 
analogue to (|3.2j) . for an integer k > do, consider the optimization problem 



(3.13) 



max rj 



s.t. w\ d = y-r}£, M k (w)h0. 
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Similarly, we choose £ to be a tms in Md(M n ) with Md (t;) y 0. Clearly, if r\ is 
feasible in (|3.13|1 . then 

(3.14) r) < V * := X min (m^)" 1 / 2 M do (y)M do ^r 1/2 ) ■ 

When rj < rf , the moment matrix Md (y — rjS,) is positive definite, and the tms 
y — rj^ is extendable to a tms w € ^ n ,2k with Mk(w) y 0. So the optimal value of 
(|3.13l) is equal to rf for every A;. Clearly, we have rf > if and only if the matrix 
Md (y) is positive semidefinite, and rf — if and only if A TO i n (-M<f (y)) = 0. 

Theorem 3.9. Let y be a tms in and £ be a tms in i?d(IR n ) with Md (£) >- 0. 

Suppose (?y, w) is feasible for \3.13]) . Let y := r/ — and fc > do. 

(i) -f/f? > and tu|t(i > d) *s /?a£, iften y belongs to Md{^ n ). 

(ii) If (i is a measure in meas(y,R n ), then there exists p € £ n ,2d smc/i i/ioi 
supp(n) CZ(p). 

In the following, suppose ji,p satisfy (ii) and {rj,w) is optimal for \3.13\ ). 

(iii) // |-Z(p)| < oo, then w\2k-2 is flat for k sufficiently large. 

(iv) If \supp(fi)\ = +oo and rankMk(w) is maximum, then w\2t can not be flat 
for allO < t < k. 

Proof, (i) is implied by y — + w\d and Theorem ll.il 

(ii) The dual optimization problem of (|3.13|) is 

(3-15) nun (/, y) s.t. (/,£} = 1, / e S„, 2do . 

Since (|3 . 13[) has a feasible w with Mk(w) y 0, the optimization problem (j3. 15[) has 
a minimizer p £ ^ n ,2d Q i and its optimal value equals rf . So 

/ pd[i= (p,y) = (p,y) -rf = 0. 

The nonnegativity of p in W 1 implies supp(//) C Z(p). 

(iii) It can be proved in a way similar to items (iii)-(iv) of Theorem 13.61 Write 
p = p\ + ■ ■ ■ +p 2 L . Then supp(/i) C Z(p) implies 

^ZpI Md {y)Pi = (p,y) = / pdfi = 0. 

Let {hi, . . . , h r } be a Grobner basis for the vanishing ideal I(Z(p)), under a total 
degree ordering. Similarly, by Theorem 12.51 for each hi , there exist I > 1 and 
/l) • • • i fr € K[x] satisfying ft? +/iPi+- • -+f r p r = 0. As in the proof of Theorem l3~6l 
we can similarly prove hi € ker Mk(w) and (fjpj, ui) = for j = 1, . . . , r, when fc is 
sufficiently large. When .2<(p) is finite, we can similarly prove w|2fc-2 is flat as for 
(iii) of Theorem EH 

(iv) The proof is same as for (iv) of Theorem 13.61 □ 

Example 3.10. Consider the same tms y as in Example 13.81 Take K = R 2 ; thus 
d g = 1 and do = 2. The moment matrix Mi(y) has rank 3 and A^iy) has rank 4, so 
y itself is not flat. The 6x6 moment matrix M2{y) is singular, so n* — 0. Choose 
£ in (|3.13p to be the tms in ^#2,4 generated by the standard Gaussian measure 
with mass one. For k — A, solve the optimization problem (|3.13|) . The computed 
optimal tms w* € .-#2,8 is not fiat, but its truncation w*\q is flat and admits a 
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unique 4-atomic measure with support {(±1,±1)}. The tms y admits the same 
measure. □ 

Unlike for the case that A is compact, the shifted tms y := y— ?y*£ does not always 
admit a measure, even when n = 1. For instance, for the tms (1, 1, 1, 1, 2) £ ^1,4, 
its moment matrix M2(y) is positive semidefmite but singular. So, 77* = but y = y 
does not admit a measure (cf. QUI Example 2.1]). 



4. A Practical Method For Finding Flat Extensions 

The preceding section discusses how to check if a tms y admits a A-measure. 
The limit Aoo of the optimal values of p. 21) or p. 51) plays a critical role. For a 
compact set A, if > 0, then y admits a A"-measure; otherwise, it does not. 
When Aoo > 0, a representing measure for y is also constructible from the relation 
y = y + Aoo£ if y admits a A-measure. In the case > 0, this approach typically 
does not give a finitely atomic representing measure, because £ is usually generated 
by a measure with infinite support. 

In many applications, one is interested in getting a finitely atomic representing 
measure. By Theorem 11.21 of Curto and Fialkow, a tms y admits a A-measure if 
and only if it is extendable to a flat tms w of higher degree. This means that if y 
admits a A-measure, then a finitely atomic measure for y can be obtained by finding 
a flat extension. In some special cases (e.g., a moment matrix is singular, or its 
associated variety is finite, etc), there exists work investigating when and how a flat 
extension could be found (cf. [7J [9J [TTJ [12] ) . For the general case, especially when 
a moment matrix is positive definite, methods for constructing a flat extension and 
determining if it exists are relatively unexplored. 

This section presents a practical SDP method for this purpose, alluded to in 
i jl.2f 2). Our numerical experiments show that it often finds a flat extension of y 
when it exists. 



4.1. TKMP. Suppose a tms w g ^ n .2k {k > d/2) is an extension of y. Clearly, 
if a truncation w\±(t > d) is flat, then the finitely atomic measure admitted by w\t 
is a representing measure for y. To find such an extension, for an integer k > d/2, 
consider the semidefmite optimization problem: 



(4.1) 



min cTw 



s.t. w\ d = y, M k -i(p* w) h 0, 
M k - dt \sf i *w) h 0, i = 0,l,. 



In the above, c k is a generic vector, and the polynomial p(x) = R 2 — \\x\W is such 
that A C B(0, R). Recall that g = 1 and M k (w) = M k ^ do (g * w). 

(1) If w* is an optimizer of (|4.1j) and w*\t(t > d) is flat, then the finitely atomic 
measure admitted by w*\t also represents y. 

(2) If y does not admit a A-measure, then (|4.ip will not be feasible for k 
sufficiently large (see the remark after Theorem 13. 3p . 

We would like to remark that the optimization problem (|4.1[) . as well as (|3.12[) , is 
an obvious and natural way to check Theorem 1 1.21 of Curto and Fialkow. 
We show some examples of illustrating this method. 
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n 


d 


d 


d 


d 


d 


2 


4(100) 


6(100) 


8(50) 


10(50) 


12(50) 


3 


4(100) 


5(100) 


6(50) 


7(50) 


8(50) 


4 


3(100) 


4(100) 


5(50) 


6(20) 


7(10) 



Table 1 . Inside the parentheses are the numbers of randomly gen- 
erated instances we tested for each pair (n, d). 



Example 4.1. Let K = jo; e : Ylt=i Xi — ^1 ■ Consider the tms y € ^#4^: 
(1; 1/5, 1/5, 1/5, 1/5; 1/10, 1/20, 0, 0, 1/10, 1/20, 0, 1/10, 1/20, 1/10). 

Its moment matrix Mi(y) is positive definite. Choose Ck to be the vector of all ones. 
Solve (|4.1j) for k = 4, and get an optimizer w* . The tms w* and its truncations 
tu*|4,u;*|6 are all flat. They admit a 5-atomic measure supported on the points: 

(0.5,0,0,0), (0,0,0,0.5), (0,0.5,0.5,0), (0,0,0.5,0.5), (0.5,0.5,0,0). 

Since w* is an extension of y, the tms y admits the same 5-atomic measure. □ 

Next we describe an experimental test of using (|4.ip to find flat extensions for 
randomly generated tms'. It produced a flat extension for all the instances we 
tested. 

Example 4.2 (random problems). Let K = [—1,1]™ be the unit hypercube. We 
generate testing instances as follows. In (|4.ip . choose Ck to be a Gaussian random 
vector, and (n, d) from Table [TJ By the theorem of Bayer and Teichmann [TJ, 
a tms y admits a if-measure if and only if y admits a iV-atomic -ftf-measure with 
N := So, we randomly generate N points u%, . . . ,um from [—1, 1]", and set 

y = a\[u\\d + • • • ajv[ujv]d for random positive numbers otj > 0. In each case, we 
solve (|4.ip starting with k = d. If an optimal w* of (|4.1j) has a flat truncation 
w*\ t {t > d), we stop; otherwise, increase k by one, and then solve (|4.ip again. 
Repeat this process until we get a flat extension of y. For each (n, d), the number 
of randomly generated instances is listed in the parenthesis after d in Table [TJ In all 
the tested instances, we got a flat truncation w*\t, for some k,t>d. Furthermore, 
the supports of the finitely atomic measures we obtained often have cardinalities 
smaller than N (we are not sure whether they are minimum or not). □ 

Our conclusion is that solving (|4.1I) is a practical method for finding a flat ex- 
tension of a tms y when it admits a ii'-measure. 

4.2. A theoretical analysis. While we do not understand well why solving (|4.1[) 
always produced a flat extension in our numerical experiments. Here, we present a 
bit of theoretical analysis suggesting that this will often produce a flat extension. 
Our analysis is to set up a linear-convex optimization problem for which (|4.1[) is an 
approximation. Then, under some assumptions, we prove that (|4.1[) asymptotically 
produces a flat extension. We start with some definitions. Denote 

Ek{y) '■= {w G .^n : 2k ■ w is feasible for (|4.1[) }. 

Eoo{y) ■= {ze^ooiK) : z\ d = y}. 

Embed Ek(y) into ^#n i0 o by the mapping w 1— > (w, 0, . . .) of adding zeros. Thus, 
Ek(y) and E oc (y) can be thought of as convex subsets of Jt n ,oo- For every w € 
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define ||u>||2 as \\w\\l — J2atN™ w a- Define the Hilbert space 

(4.2) Jil^ := {w e J? n ,oc ■ \\w\U < oo}- 

Up to a shifting and scaling, we can generally assume K C B(0,R) with R < 1. 

We assume the vectors in (|4.1I) for different k are consistent, i.e., each Ck-i is 
a leading subvector of c/.. Write Ck(x) — c£[x]2k- Then, its limit is a real analytic 
function, which is denoted as c(x) = T, a c a x a . The analytic function c(x) can also 
be thought of as a vector c G s^n.oa- If || C I|2 < o°j then c is a continuous linear 
functional acting on ^ as 

(c,w) := y^c a w a . 

a 

Now we consider the optimization problem: 

(4.3) min / c(x)dp, s.t. fi € meas(y, K). 

P Jk 

When y belongs to Md{K), its feasible set is nonempty. The objective is a linear 
functional acting on the convex set meas(y, K). Note that a vector w belongs to 
Eoaiy) if and only if it admits a measure in meas(y, K). Thus, the optimization 
problem (|4.3[) is also equivalent to 

(4.4) min (c, w) s.t. w € E oc (y). 

w 

If the optimization problem (|4.3[) has a unique minimizer /i*, then /i* must have 
finite support (see Appendix). Thcorcm l4. 31 below shows that if the biggest support 
of minimizing measures of (|4.3j) is finite, then (|4.1[) asymptotically yields a flat 
extension as k goes to infinity. Theorem A. 2 in Appendix shows that for a dense 
subset of analytic functions c in ^{K ) (the space of all continuous functions defined 
on K), this condition (see (|4.5|) ) holds. 



Theorem 4.3. Assume K C B(0,R) with R < 1, the c' k s are consistent, and c is 
in ^#^ oc . Lei Sol be the set of all optimizers of f^. 3\ ). Suppose there exists N > 
such that 

(4.5) j (J {supp([i*)}\<N. 

Let w/ fc ) &e an optimizer of j4-l\ l. Then, for every r > d g N (d g is given by U.9\) ) 
and every /i* € Sol, the tms z* := J K [x]2 r dfJ,* is fiat and 

(4.6) lim dist(w^\ 2r , U ) = w/iere U=\ [x] 2r dfi* : ^* S Sol 

In particular, if Sol is a singleton, then {w' fe '|2r} converges to aflat tms. 

Proof. Let ro := d g iV. If Zq = 0, then z* must be a zero vector, and it is clearly 
flat. So, we consider the general case Zq > 0. We claim that a truncation z*\i of z* 
must be flat for some < I < ro- Otherwise, if not, then we must have 

1 = rankAf (z*) < rank M dg (z*) < rank M 2dg (z*) < ■ ■ ■ < rankM r „(z*). 

The above then implies the contradiction: 

|supp(//)| > rankAf ro (z*) > 1 + r /d g > N. 
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So, z*\i is flat for some I < r$. By Theorem we know fi* is the unique rep- 
resenting measure of z*\i. Since z* is an extension of z*\i and admits the same 
measure /«*, we know that z* must also be flat. 

We prove (|4.6p by contradiction. Suppose otherwise it is false, then 

(4.7) dist(w (fc) | 2 r, U) > a for some a > 0. 

By Lemma rA.31 we know the sequence {w^ } is bounded in the Hilbert space ao . 
By the Alaoglu Theorem (cf. [2] Theorem V.3.1]), {w^} has a subsequence that is 
convergent in the weak-* topology. That is, it has a subsequence, which we denote 
again by {w^} for convenience, such that for some w* € j^^, 

£{w {k) ) ^ l{w*) 

for every continuous linear functional in ^ (the space x is self-dual) . This 
implies (e.g., choose i as l(w) — w a for each a <E N n ) that for all t 

(4.8) w {k) \ 2t -> w*\ 2t , and w*\ d = y. 
Since u;W is feasible in (14.11) . from the above we get that 

M t -d i (gi*w*) t 0(i = 0,l,...,m), M t _i(p*io*) h 0. 

Hence, the vector u>* is a full moment sequence whose localizing matrices of all 
orders are positive semidefmite (because we have w*\ 2t € Qk(K)* for all t). By 
Lemma 3.2 of Putinar [23], we know w* admits a i^-measure and w* belongs to 
Eoo(y)- Note that c is a continuous linear functional acting on ^ as (c, u>). So, 

(c,iu*) = lim (c, w; (fc) ). 

By Lemma [A. 41 we know that is also an optimizer of (|4.4[) and w*|2r belongs to 
J7. However, (|4.8j) implies the convergence w/^^r — > w*|2r, and (|4.7|) implies that 
w*\ 2r does not belong to U, a contradiction. Therefore, (14.61) must be true. 

When Sol is a singleton, {w^\ 2r } clearly converges to a flat tms. □ 

By Theorem l4.3[ we know that every limit point of the sequence {w^ \ 2r } is flat. 
We would like to remark that the conditions in Theorem |L3] are only for theoretical 
analysis. In practical implementations, they are not required. 

4.3. The case K = R". When K = W 1 , we propose a similar method like (I4.1[) 
for solving TMPs. To find a finitely atomic representing measure for a tms, it is 
enough to find one of its flat extensions. Like (14. ip . a practical method is that for 
an order k > d/2 we solve the SDP problem: 

{min cTw 
s.t. w\ d = y, M k {w) h 0. 

For k > d/2, the above feasible set is typically unbounded. To guarantee (|4.9j) has 
a minimizer, we usually choose c k as c[[x] 2 k — [x]JCk[x]k with Ck h 0. A simple 
choice for Ck is the identity matrix. 

Example 4.4. Consider the tms y £ ^#2,4: 

(1; -1/6, 1/2; 3/2, 0, 3/2; -7/6, 1, 1/3, 3/2; 7/2, -1, 2, 1, 7/2). 

Its moment matrix M 2 {y) is positive definite. By Theorem 3.3 of |13) . this tms y 
admits a measure. Choose Ck as c^[x] 2 k — [x]^[x]k, and solve (|4.9|) for k = 4. The 
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computed optimal w* 
the points 0: 



2,8 is flat. It admits a 7-atomic measure supported on 

0.4613), 



□ 



(-1.1902, -1.4545), (-0.2725, -1.4977), (1.4406, -0.9396), (0.2150 
(-1.7147, 0.1952), (-1.6942, 1.1593), (0.9645, 1.7098). 

The tms y admits the same 7-atomic measure. 

We would like to remark that solving (|4.ip or (|4.9p might not give a representing 
measure for y whose support has minimum cardinality. For instance, in Exam- 
ple ^. 4| the tms y there admits a 6-atomic measure (with support {±(1, 1), ±(1, —1), 
(1, 2), (—2, 1)}, but we got a 7-atomic measure by solving (|4.9|) . 

Example 4.5 (random examples). Now we test the performance of solving (I4.9[) 
for finding flat extensions. We choose Ck such that cjw — Ck • Mk(w) with Ck = 
LL T and L being a Gaussian random matrix. Select (n, d) from Table [TJ Like in 
Example 14. 2 [ for each pair (n, d), randomly generate N = ( n ~ d ~ d ) points u\,...,un 
obeying a Gaussian distribution, and set y = ai[wi]d + ■ ■ ■ aj^[uiq] d (at > is 
randomly chosen). For each instance, solve (|4.9[) starting with k — d. If an optimal 
w* of (|4.9[) has a flat truncation w*\ t {t > d), we stop; otherwise, increase k by 
one, and solve (|4.9I) again. Repeat this process until a flat extension of y is found. 
In all the tested instances, we were able to find a flat extension of y for some 
k > d. In the computations, we sometimes obtained |supp(/x)| < ./V and sometimes 
|supp(/i)| > N, while the former was gotten more often. □ 

The above random example leads us to believe that (|4.9[) is practical for solving 
TMPs. We now give a sufficient condition for (|4.9I) to produce a flat extension. 

Theorem 4.6. In let d = 2do + 1, k = do + 1 and c[[x]2k = [x]^ Ck[x]k be 

such that Ck y 0. If the tms y 6 *$ n ,d admits a rank Md (y)- atomic measure, then 
the optimizer w* of J^.ff| ) is flat. 

Proof. Let r = rank Md (y) and fi be a r-atomic measure admitted by y 



Set 



z = [x] d+ idn, M do+1 (z) 



M do (y) 
B T 



C k 



C n 
c T 

°12 



C12 
C22 



The matrix B consists of moments y a with |a| < d. Note that 

r = |supp(/x)| > rankA/ do+ i(z) > rankM do (y) = r. 

So rank Md„+i(z) = r. Then Md„+i(z) y implies 

Z = B T M do (y) + B>{), 

where the superscript + denotes the Moore-Penrose Pseudo inverse of a matrix. 
For every w that is feasible for (|4.9p . write 

~M do (y) B 
B T W 



A! 



d + l 



C k • M do+1 {w) = C n • M do (y) + 2C 12 •B + C 22 »W. 
Thus, the objective of (14.9[) can be equivalently replaced by C 2 2*W . If M do +i(w) y 
0, then W - Z y and C 22 • W > C 22 • Z. So, C 22 • Z is a lower bound of C 22 • W 
for all feasible w. This implies that z is an optimizer of (|4.9[) . Indeed, z is the 
unique optimizer of (|4.9[) . Suppose w is another optimizer, then C 22 • (W — Z) = 0. 



^Here only four decimal digits are shown to present points of the support. 
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mm c k w 



The fact that W > Z and C22 >~ implies W = Z. Hence w — z. Clearly, z is flat, 
and the proof is complete. □ 

Remark: When K in (jl.ip is noncompact, we can similarly solve the SDP problem: 
(4-10) 

{ s.t. w\ d = y,M t ^ di (gi*w)hO,i = 0,l,...,m 

in analogue to (|4.1j) and (|4.9|) . To guarantee it always has a minimizer w*, we 
typically choose as c[[x]2k = Ck[x]k with Ck >~ 0, like in (|4.9|) . If a truncation 
u>*| t (t > d) is flat, we get a flat extension of y. 

Acknowledgement The authors thank Larry Fialkow for helpful comments on 
this work. J. William Helton was partially supported by NSF grants DMS-0700758, 
DMS-0757212, and the Ford Motor Co. Jiawang Nie was partially supported by 
NSF grants DMS-0757212 and DMS-0844775. 

Appendix A. Optimization with measures 



We prove that for a dense set of optimization problems (J473J) , the optimal solu- 
tions are measures having finite supports. This justifies condition (|4.5|) in <j47 



Lemma A.l. Let y be a tms in J% n ,d and K be a subset o/K™. Assume the set 
meas(y, K) is nonempty. If a measure [i is an extreme point of meas(y, K), then 
fi is finitely atomic and | supp(fi) \ < (™ J d ) • 

Proof. We prove by contradiction. Suppose n is extreme and |supp(/i)| > N := 
(™ J d ) ■ Then we can decompose supp(/i) in a way such that 

JV+l 

supp(M) = |J Si, S l nS J =9, fi(S j )>0, Wi^j. 
J=l 

For each j, let jij be the restriction of n on Sj , i.e., fij = . Note that dim R[x]d = 
N. So there exists (ai, . . . , ajv+i) 7^ satisfying 



N+l 

E 



a,j / x a d^.j = 0, V|a| < d. 
Jk 



ik 

Define a measure v (possibly non-positive) such that 

v = u\si H h ^|s«+i, Hs* = ctjfij, j — 1, . . . , N + 1. 

Note that supp(^) C supp(ju). Then, for e > sufficiently small, /i ± ev are 
nonncgative Borel measures supported on K , representing the same tms y and 

H = (n + eu)/2+(jj,-ev)/2. 

This contradicts the extremality of /1 in meets (y, K). □ 

Theorem A. 2. Let y be a tms in ^ n ,d- For any real analytic function f on a 
compact set K and e > 0, there exists an real analytic function f on K such that 
11/ — /II 00 < £ an d every minimizer of 

(A.l) min / fdfi s.t. fi G meas(y, K) 

t* Jk 
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has cardinality at most • Let Sol be the set of all optimizers above. Then 
(A.2) | |J {su PP (n)}\< f n+ d d ). 

Proof. For a given real analytic function / on K, consider the optimization problem 

(A.3) min / fd[i s.t. fi G meas(y, K). 

* Jk 

The set of positive measures on K whose total masses equal yo is compact in the 
weak-* topology (denote this topology by T) by the Alaoglu Theorem (cf. [2, 
Theorem V.3.1]). Recall that the weak-* topology is the topology on the measures 
regarded as the dual space of the Banach space ^(K). It is the weakest topology 
for which the convergence of a sequence of measures /j, n — > /j, implies that for every 



(A.4) / hdfi n — !• / hdfi. 

Jk Jk 

This implies that every moment of fj, n converges to the corresponding one of /i. 
Hence, meas(y, K) is 7~-closed inside a compact set, and it is also T-compact. 

Thus, by compactness of its feasible set, the optimization problem (|A.3j) has a 
minimizer, which is a generalized convex combination of certain extreme points, by 
the Choquet-Bishop-de Leeuw Theorem (cf. [33]). To be more precise, this means 
that for every fi € measijj, K), there exists a probability measure T on the set £ of 
extreme points of meas(y, K) such that 

l(fi)=J i(ji) dT(/x) 

for every affine function £ on meas(y,K). 

Let 7* be the optimal value of (|A.3[) and fx be a minimizer, and let T denote the 
probability measure on £ representing jl and let £ denote the support of Y. Then 



By optimality of 7*, /d/i > 7* for all fi E £. Indeed, J K fdfi — 7* for all /j, G £ . 
Otherwise, suppose J K fdfi > 7* on a set of /j, having positive T measure. Then 

7* = J(j K fd f ?\df( f i)>j*, 

which yields a contradiction. This implies J K fdfi = 7* on £. Choose a measure 
/i* € £. It is extreme and the optimum of (|A.3I) is attained at n* . 

By Lemma fA. 11 we have |supp(/i*)| < N :— ( n ^ d ). Without loss of generality, 
we can normalize yo = J K d/i* = 1 and denote supp(/i*) = {ui, . . . , u r }. Let e(x) 
be the exponential function defined as 

e(x) = e • exp(— ||x — U1II2 •• ■ — ttr- 1 1 1) ■ 

Clearly, e(x) = e for all x € supp(/i*) and < e(x) < e for all x $ supp(^*). This 
implies that 



max / e c?/i = / e dp,* 

IJ.emeas(y,K) J K J K 
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Set / := / — e and 

(A. 5) (3 := min / / dfi. 

fiemeas{y,K) J K 

Then 

(A. 6) f3 > min / / dfi — max / e dfi = / / d/x* — e. 

fj,emeas(y,K) J K ^emeas(y,K) J K J K 

On the other hand 

/ dfi* — e = f dfi* > min / f dfi = (3 

K JK n£meas(y,K) J K 

because fj,* belongs to meas(y, K). Thus 

P= I fdfx*-e = [ fdfx* 

JK JK 

and n* is a minimizer of (|A.5|) as well as (|A.1[) . 

Suppose jl* is another minimizer of (|A.5|) . We wish to show that 

(A. 7) suppt(jl*) C {u\, . . . , u r }. 

If this is not true, then J e dfl* < e and 



(3= f dfi* — e < min / / dfi — e dp,* . 
JK ii£meas(y,K) J K J K 

This implies 

<3< / f dfi* - e dfi* = / fdfl* = i3, 

JK JK JK 

which is a contradiction. 

The inequality (|A.2[) follows from the first part, because the average of all the 
minimizers is still a minimizer. □ 

Remark: As we can see in the above proof, if fi* is a unique minimizer of (|4.3I) . 
then fi* must have finite support and |supp(^*)| < ( n ^ d )- 

The following lemmas are used in the proof of Theorem 14.31 

Lemma A. 3. Suppose K C _B(0, R) with R < 1 and c is in ^? loo . If a tms w 
belongs to Ek(y) U E^y), then 

lk|| 2 <y /(i-i? 2 ), 

and for any integer t > it holds that 

(A.8) | ]T c a w a \ < \\c\\ 2 ■ y ■ R 2t+2 /(l - R 2 ). 

\a\>2t 

Proof. First, we consider the case that w € Ek(y). The condition M t -i(p * w) >z 
in (|4.1|) implies that for every t = 1, . . . , k 

r 2 j?U\\x\\ 2 2~ 2 )-^U\\x\\¥)>o- 

A repeated application of the above gives 

^(IMI^) < R 2t yo- 
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Since Mk(w) is positive semidefmite, we can see that 

fe 

IMIa < Wk{w)\\ F < Trace(M k (wj) = E E W2a > 

t=Q \a\=t 

E w ^ = E x2a ) ^ ^(Mll*) < R 2t vo- 

\a\=t \a\=t 

Thus, it holds that 

1Mb < y (i + R 2 + ■ ■ ■ + R 2k ) < vol {i - R 2 )- 

When k > t, we have 

2 



uu 2 , 



E c «H ^( E c ')( E 

2t<\a\<2k 2t<\a\<2k 2t<|a|<2fc 

Let M t ^k be the principle submatrix of Mk(w) whose row and column indices (3 
satisfy t < \/3\ < k. Clearly, M t ^ cl and we similarly have 

1 /2 

( E <\Wt,k\\F<Trace{M t , k ) 

2t<\a\<2k 

k k k 

= E E ^ a < E ^m\1 3 ) < E y» R2j - 

j=t+l \ a \-j j=t+l 3 =t+l 

Combining all of the above , we get 

E c a w a \<\\c\\ 2 - yo -(R 2t+2 + --- + R 2k )<\\c\\ 2 -yo-R 2t+2 /(l-R 2 ). 

2t<\a\<2k 

Note that w a = if \a\ > 2k. So, (|A.8|) is true. 

Second, we consider the case that w £ E oa (y). The sequence w admits a K- 
measure /i such that w a = J x a d[i for every a. For every k, consider the truncation 
z = J K [x]2kdfJ,. Clearly, the tms z is feasible for (|4.ip . By part (i), we have 

\\z\\ 2 <y /(l-R 2 ), \ E c a z a \<\\c\\ 2 -y Q -R 2t+2 /(l-R 2 ). 

2t<\a\<2k 

The above is true for all k, and implies (|A.8[) by letting k —> oo. □ 

Lemma A. 4. Assume K C 5(0, i?) u>i£/i i? < 1, and c is m oc . Let 7^,7 &e 
i/ie optimal values of |^.i| ) and respectively. Then, we have the convergence 

7^ — > 7 as — > 00. 



Proof. By Lemma I A. 3 [ for an arbitrary e > 0, there exists t such that 
(A.9) J ^ c a w Q < e, VttieBk(y)U E^y). 

\a\>t 

This implies that for every w £ E k {y) U E oc (y) with > t, it holds that 

(A.10) \(c,w)-cjw\ 2t \<e. 

Now we consider the truncated optimization problems 

(A. 11) min cjw\ 2t s.t. w £ E k (y) 
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and 

(A.12) min cFwkt s.t. w € E^y). 

w 

Let Tfc and r be the optimal values of (|A. 1 1|) and (|A.12|) respectively. Note that 
(|A.11| is a semidefinite relaxation of (|A.12|) . Thus, we have — > r, by Theorem 1 
of Lasserre [T7] . The inequality (|A.10I) implies that for all k > t 

min cjw\2t — e < min w < min cjw\2t + e- 

wGE k (y) w£E k (y) weE k (y) 

This shows that \jk — t^| < e for k > t. Applying a similar argument to (|A.12|) . we 
can get |7 — t| < e. Note that 

\lk - 7l < l7fc - T fc| + l r fe - T l + l r - 7l- 
Because Tk r, one can get that 

< limsup \ jk — 7| < 2e. 

Since e > is arbitrary, we must have 7fc — > 7 as fc — > 00. □ 
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